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ABSTRACT 

Context. The first detection of a moon around an extrasolar planet (an “exomoon”) might be feasible with NASA’s Kepler or ESA’s 
upcoming PLATO space telescopes or with the future ground-based European Extremely Large Telescope. To guide observers and to 
use observational resources most efficiently, we need to know where the largest, most easily detected moons can form. 

Aims. We explore the possibility of large exomoons by following the movement of water (H 2 0) ice lines in the accretion disks around 
young super-Jovian planets. We want to know how the different heating sources in those disks affect the location of the H 2 0 ice lines 
as a function of stellar and planetary distance. 

Methods. We simulate 2D rotationally symmetric accretion disks in hydrostatic equilibrium around super-Jovian exoplanets. The 
energy terms in our semi-analytical framework - (1) viscous heating, (2) planetary illumination, (3) accretional heating of the disk, 
and (4) stellar illumination - are fed by precomputed planet evolution models. We consider accreting planets with final masses between 
1 and 12 Jupiter masses at distances between 1 and 20 AU to a solar type star. 

Results. Accretion disks around Jupiter-mass planets closer than about 4.5 AU to Sun-like stars do not feature H 2 0 ice lines, whereas 
the most massive super-Jovians can form icy satellites as close as 3 AU to Sun-like stars. We derive an empirical formula for the 
total moon mass as a function of planetary mass and stellar distance and predict that super-Jovian planets forming beyond about 
5 AU can host Mars-mass moons. Planetary illumination is the major heat source in the final stages of accretion around Jupiter-mass 
planets, whereas disks around the most massive super-Jovians are similarly heated by planetary illumination and viscous heating. This 
indicates a transition towards circumstellar accretion disks, where viscous heating dominates in the stellar vicinity. We also study a 
broad range of circumplanetary disk parameters for planets at 5.2 AU and find that the H 2 0 ice lines are universally between about 15 
and 30 Jupiter radii in the final stages of accretion when the last generation of moons is supposed to form. 

Conclusions. If the abundant population of super-Jovian planets around 1 AU formed in situ, then these planets should lack the 
previously predicted population of giant icy moons, because those planets’ disks did not host H 2 0 ice in the final stages of accretion. 
But in the more likely case that these planets migrated to their current locations from beyond about 3 to 4.5 AU they might be orbited 
by large, water-rich moons. In this case, Mars-mass ocean moons might be common in the stellar habitable zones. Future exomoon 
detections and non-detections can provide powerful constraints on the formation and migration history of giant exoplanets. 

Key words. Accretion, accretion disks - Planets and satellites: formation - Planets and satellites: gaseous planets - Planets and 
satellites: physical evolution - Astrobiology 


1. Introduction 


Now that the detection of sub-Earth-sized objects has become 
possible with space-based photometry (Muirhead et al. 2012; 
Barclay et al. 2013) and almost 2000 extrasolar planets have 
been confirmed (Batalha et al. 2013; Rowe et al. 2014), tech¬ 
nological and theoretical advances seem mature enough to find 
moons orbiting exoplanets. Natural satellites similar in size to 
Mars (0.53 Earth radii, R®) or Ganymede (0.41 R®) could be 
detectable in the available Kepler data (Kipping et al. 2012; 
Heller 2014), with the upcoming PLATO space mission (Simon 
et al. 2012; Rauer et al. 2014) or with the European Extremely 
Large Telescope ( E-ELT ) (Heller & Albrecht 2014). Exomoons 
at about 1 AU or closer to their star might be detected during stel¬ 
lar transits (Sartoretti & Schneider 1999; Szabo et al. 2006; Kip¬ 
ping et al. 2009). Young, self-luminous planets beyond 10 AU 
might reveal their satellites in the infrared through planetary 
transits of their moons, if the unresolved planet-moon binary 
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can be directly imaged (Peters & Turner 2013; Heller & Al¬ 
brecht 2014). Moon formation theories can guide observers and 
data analysts in their searches to streamline efforts. In turn, the 
first non-detections of exomoons (Brown et al. 2001; Pont et al. 
2007; Kipping et al. 2013b,a, 2014) and possible future findings 
provide the first extrasolar observational constraints on satellite 
formation. 

In a recent study, we developed a circumplanetary disk accre¬ 
tion model and simulated the radial motion of the water (H 2 0) 
ice line around super-Jovian 1 exoplanets at 5.2 AU around Sun¬ 
like stars (Heller & Pudritz 2014). The H 2 O ice line is critical 
for the formation of large, possibly detectable moons, because 
here the surface density of solids increases sharply by a factor of 
3 to 4 (Hayashi 1981). The composition and the masses of the 
Galilean moons are usually considered as records of the location 
of the H 2 O ice line and, more generally, of the temperature distri- 


1 Throughout the paper, our reference to “super-Jovian” planets in¬ 
cludes planets with masses between that of Jupiter (M Jllp ) and 12 M Jup , 
where the latter demarcates the transition into the brown dwarf regime. 


Article number, page 1 of 11 





A&A proofs: manuscript no. ms 


bution in Jupiter’s accretion disk (Pollack & Reynolds 1974; Lu- 
nine & Stevenson 1982). The inner moons Io and Europa turned 
out mostly rocky and comparatively light, whereas Ganymede 
and Callisto formed extremely rich in water ices (about 50 %) 
but became substantially more massive (Showman & Malhotra 
1999). 

Our main findings were that (1) super-Jovian planets at 
5.2 AU around Sun-like stars have their JOO ice lines between 
about 15 and 30 Jupiter radii (Rj up ) during the late stages of ac¬ 
cretion, when the final generation of moons form. This range 
is almost independent of the planetary mass (M p ). (2) With the 
most massive planets having the most widely extended disks, 
these disks host the largest reservoirs of JTO ices. In particular, 
the total instantaneous mass of solids (M s y) in these disks scales 
roughly proportional with M p . (3) The current orbital position 
of Ganymede is very close the mean radial location of the cir- 
cumjovian JTO ice line in our simulations, suggesting a novel 
picture in which Ganymede formed at a circumplanetary ice line 
trap. (4) Heat transitions, however, transverse the accretion disks 
on a very short timescale (~ 10 4 yr) once the planet opens up 
a gap in the circumstellar accretion disk and thereby drastically 
reduces the supply of material. This timescale is short compared 
to the satellite migration time scale (10 5 - 10 7 yr) (Tanaka et al. 
2002; Canup & Ward 2002, 2006; Sasaki et al. 2010; Ogihara 
& Ida 2012) and the time that is left for the final accretion until 
shutdown (10 5 yr). Hence, heat transitions around super-Jovian 
planets cannot act as moons traps, which indicates a different be¬ 
havior of planet and moon formation (Menou & Goodman 2004; 
Hasegawa & Pudritz 2011; Kretke & Lin 2012). 

We here extend the parameter space of our previous paper 
and consider planetary accretion between 1 and 20 AU from the 
star. This range is motivated by two facts. First, super-Jovians are 
the most abundant type of confirmed planets, that is, objects with 
accurate mass estimates at these distances (see Fig. I) 2 . And sec¬ 
ond, this range contains the stellar habitable zone around Sun¬ 
like stars, which is of particular interest given the fact that gi¬ 
ant, water-rich moons might be frequent habitable environments 
(Heller & Armstrong 2014; Heller et al. 2014). Our main goal 
is to locate the H 2 O ice lines in the circumplanetary accretion 
disks at the time of moon formation shutdown. Their radial sep¬ 
arations from the planet will correspond to the orbital radii at 
or beyond which we can suspect the most massive, water-rich 
moons to reside. In particular, stellar illumination at only a few 
AU from the star will prevent the accretion disks from having 
H 2 O ice lines and therefore from hosting large moons. We deter¬ 
mine these critical stellar distances below. 


2. Methods 

We use the framework developed in Heller & Pudritz (2014), 
which models a 2D axisymmetric accretion disk in hydrostatic 
equilibrium around the planet. It considers four heating terms, 
or energy fluxes, of the disk, namely (1) viscous heating (/-\.j s ), 
(2) planetary illumination ( F p ), (3) accretional heating of the 
disk (F acc ), and (4) stellar illumination, all of which determine 
the midplane and surface temperature of the disk as a function 
of planetary distance (r). This model (based on earlier work 
by Makalkin & Dorofeeva 1995, 2014) considers F v j s as a dis¬ 
tributed energy source within the disk, while F p , F acc , and stel¬ 
lar illumination are considered external heat sources. The model 

2 Data from the Kepler space telescope suggests, however, that terres¬ 
trial planets are more abundant at about 1 AU around Sun-like stars than 
gas giant planets (Dong & Zhu 2013; Mulders et al. 2015). 



Fig. 1. Stellar distances and planetary masses of extrasolar planets listed 
on www. exoplanet. eu as of 5 April 2015. Symbols indicate the dis¬ 
covery method of each planet, six solar system planets are shown for 
comparison. Note the cluster of red dots around 1 AU along the abscissa 
and between 1 and 10 Mj up along the ordinate. These super-Jovian plan¬ 
ets might be hosts of Mars-mass moons. The shaded region denotes the 
solar habitable zone defined by the runaway and maximum greenhouse 
as per Kopparapu et al. (2013). 


also includes an analytical treatment for the vertical radiative en¬ 
ergy transfer, which depends on the Planck mean opacity (/fp). In 
real disks, /fp will depend on the disk temperature and the com¬ 
position of the solids. In particular, it will be a function of the 
radial distance to the planet and it will evolve in time as small 
particles stick together and coagulate. To reduce computational 
demands, we here assume a constant /fp throughout the disk, but 
we will test values over two orders of magnitude to explore the 
effects of changing opacities. The vertical gas density in the disk 
is approximated with an isothermal profile, which is appropri¬ 
ate because we are mostly interested in the very final stages of 
accretion when the disk midplane and the disk surface have simi¬ 
lar temperatures. Furthermore, deviations between an isothermal 
and an adiabatic vertical treatment are significant only in the very 
dense and hot parts of the disk inside about 10Rj up . As the H 2 O 
ice line will always be beyond these distances, inaccuracies in 
our results arising from an isothermal vertical model are negligi¬ 
ble. 

The heating terms (l)-(3) are derived based on precomputed 
planet evolution models (provided by courtesy of C. Mordasini, 
Mordasini 2013), which give us the planetary mass, planetary 
mass accretion rate ( M ), and planetary luminosity ( L p ) as a func¬ 
tion of time (f) (see Fig. 1 in Heller & Pudritz 2014). L p (t) is a 
key input parameter to our model as it determines F p (t), and it 
is sensitive to the planet’s core mass. Yet, in the final stages of 
planetary accretion, L p differs by less than a factor of two for 
planetary cores between 22 to 130 Earth masses. With F p oc L p 
and temperature scaling roughly with (Fp) 1 ^ 4 , uncertainties in 
the midplane or surface temperatures of the disk are lower than 
20 %, and uncertainties in the radial distance of the H 2 O ice line 
are as large as a few Jovian radii at most. For all our simulations 
the precomputed planetary models assume a final core mass of 
33 Earth masses, which is about a factor of three larger than the 
mass of Jupiter’s core (Guillot et al. 1997). With higher final core 
masses meaning higher values for L p at any given accretion rate. 
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Table 1. Parameterization of the circumplanetary disk as described in Heller & Pudritz (2014). 


Symbol 

Meaning 

Fiducial Value 

Constant or parameterized planetary disk parameters 

r 

distance to the planet 


Rh 

planetary Hill radius 


r c f 

centrifugal radius 


A 

transition from optically thick to optically thin (viscously spread) outer part 

27.22/22 x r cf (fl) 

r& 

outer disk radius 

R h / 5 (i ’ e) 

A (r) 

radial scaling of S(r) 


T m (r) 

midplane temperature 


T s (r) 

surface temperature 


I(r) 

gas surface density 


h(r) 

effective half-thickness 


Po(r) 

gas density in the midplane 


Ps(r) 

gas density at the radiative surface level 


Zs(r) 

radiative surface level, or photospheric height 


q s (r) 

vertical mass coordinate at the radiative surface 


c s (T m (r)) 

speed of sound 


a 

Viscosity parameter 

0.00 \ {d) 

v(r, T m ) 

gas viscosity 


r 

adiabat exponent, or ratio of the heat capacities 

1.45 

X 

dust enrichment relative to the protostellar, cosmic value 

10 {c) 

A d 

dust-to-mass fraction 

0.006^ 

P 

mean molecular weight of the H/He gas 

2.34 kg/mol^ 

Variable planetary disk parameters 

Kp 

Planck mean opacity 

10~ 3 - 10~‘ m 2 kg^ 1 ® 

k s 

fraction of the solar radiation flux contributing to disk heating at z < z s 

O 

1 

O 

'c? 

^shut 

shutdown accretion rate for moon formation 

100,10,1 M Gan /Myr 


Notes. See Heller & Pudritz (2014) and additional references for details about the respective parameter values or ranges: (a) Machida et al. (2008) 
{b) Sasaki et al. (2010) (c) Makalkin & Dorofeeva (2014) (rf) Keith & Wardle (2014) (e) Lunine & Stevenson (1982) Hasegawa & Pudritz (2013) 
(s) D'Angelo & Bodenheimer (2013) (h) Bell et al. (1997) . 


we actually derive upper, or outer, limits for the special case of 
the HiO ice line around a Jupiter-like planet 5.2 AU around a 
Sun-like star. 

We interpolate the Mordasini (2013) tracks on a linear scale 
with a time step of 1 000 yr and evaluate the four heating terms as 
a function of r, which extends from Jupiter’s co-rotation radius 
(2.25 /(j up in our simulations, Sasaki et al. 2010) out to the disk’s 
centrifugal radius (r c f). At that distance, the centrifugal force act¬ 
ing on a gas parcel equals the gravitational pull of the planet. We 
compute r c f using the analytical expression of Machida et al. 
(2008), which they derived by fitting a power law expression to 
their 3D hydrodynamic al simulations of circumplanetary accre- 
tion disks. In this model, r c f oc for super-Jovian plan¬ 

ets at a given stellar distance. Viscous heating is governed by 
the a viscosity parameter (Shakura & Sunyaev 1973), which we 
fix to a value of 10~ 3 in our simulations (for a discussion see 
Sect. 4 in Heller & Pudritz 2014). Note that even variations of 
a by an order of magnitude would only change our results for 


the circumplanetary H 2 O ice line location during the final stages 
of planetary accretion by a few planetary radii at most, because 
then the disk is mostly heated by planetary illumination. 

The sound velocity (c s ) in the disk midplane is evaluated 
as c s = 1.9kms _1 V7m(r)/1000K (Keith & Wardle 2014), 

which is an adequate approximation since ionization can be ne¬ 
glected in the late stages of moon formation when disk temper¬ 
atures are usually below 1 000 K. At each time step, we assume 
a steady-state gas surface density (Z„) that is derived analyti¬ 
cally by solving the continuity equation of the mass inflow onto 
a centrifugally supported, viscous disk with a uniform flux per 
area (Canup & Ward 2002). The equations of energy transport 
(Makalkin & Dorofeeva 2014) then allow us to derive the tem¬ 
perature profile both at the disk surface, where the optical depth 
r = 2/3, and in the disk midplane. This model invokes ab¬ 
sorption of planetary illumination in the disk photosphere at a 
height z s (r) above the midplane, modeled by an absorption co¬ 
efficient ( fc s ), as well as the transport of energy through an opti¬ 
cally thick disk to the surface, modeled by Kp. For the dust-to-gas 
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Fig. 2. Temperature structure in the disk around a Jupiter-mass planet 5.2 AU from a Sun-like star. Solid black lines indicate disk midplane 
temperatures, solid gray lines disk surface temperatures. The other lines indicate a hypothetical disk surface temperature assuming only one heat 
source (see legend), (a): At 0.992 Myr in this particular simulation, the disk has sufficiently cooled to allow the appearance of an H 2 0 ice line in the 
disk midplane at the outer disk edge at 25.8 _ftj up . The planetary accretion rate is about 4 x 10 2 M Gan Myr -1 . (b): At 1.073 Myr, when the planetary 
accretion rate has dropped to 10 M Gan Myr -1 , the H 2 0 ice line has settles between the orbits of Europa and Ganymede (see colored symbols). In 
these final stages of accretion, disk temperatures are governed by planetary illumination (see dashed line for T p ). 





Fig. 3. Same as Fig. 2, but now for a 12Mj up planet. Note that the accretion disk is substantially larger than in the 1 Mj up case, (a): The H 2 0 ice 
line appears much later than in the Jovian scenario, here at 1.334 Myr. The planetary accretion rate is about 10 3 M Gan Myr -1 and viscous heating is 
dominant in this phase, (b): Once the planetary accretion rate has dropped to 10 M Gan Myr -1 at 1.43 Myr, the H 2 0 ice line is at a similar location 
as in the Jupiter-like scenario, that is, between the orbits of Europa and Ganymede. 


ratio (Xd) we take a fiducial value of 0.006 (Lunine & Steven¬ 
son 1982; Hasegawa & Pudritz 2013) in the inner dry regions 
of the disk, and we assume that it jumps by a factor of three at 
the H 2 O ice line. The resulting mathematical framework con¬ 
tains several implicit functions, which we solve in an iterative, 
numerical scheme. Table 1 gives a complete overview of the pa¬ 
rameters involved in the model. Our model neglects the effects 
of planetary migration, which is an adequate approximation be¬ 
cause the planets don’t migrate substantially within the w 10 5 yr 
required for moon formation. We also do not actually simulate 
the accretion and buildup of moons (Ogihara & Ida 2012). 

As an extension of our previous study, where all planets were 
considered at 5.2 AU around a Sun-like star, we here place hy¬ 
pothetical super-Jovian planets at stellar distances between 1 and 
20 AU to a solar type host. Although the precomputed planet for¬ 
mation tracks were calculated at 5.2 AU around a Sun-like star. 


we may still consider different stellar distances because the ac¬ 
cretion rates through any annulus in the circumstellar disk are 
roughly constant at any time (and thereby similar to the instan¬ 
taneous stellar accretion rate). The surface densities of gas and 
solids are naturally lower at larger stellar separations. Hence, the 
accretion rates provided by the tracks will actually overestimate 
M and L p in wider orbits at any given time. However, this will 
not affect our procedure as we are not primarily interested in the 
evolution as a function of absolute times but rather as a func¬ 
tion of accretion rates. In particular, we introduced the shutdown 
accretion (M s h u t) in Heller & Pudritz (2014), which we use as a 
measure to consider as the final stages of moon formation around 
accreting giant planets. M s h ut is a more convenient quantity (or 
independent variable) than time to refer to, because, first, the ini¬ 
tial conditions of planet formation as well as the age of the sys¬ 
tem are often poorly constrained. And, second, accretion rates 
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Fig. 4. Surface densities around the four test planets shown in Figs. 2 
and 3. Black lines refer to the Jupter-mass planets, gray lines to the 
12iWj up -mass super-Jovian. Note that the disk of the Jupiter twin is 
much smaller than that of the super-Jovian and that both surface density 
distributions decrease with time. The radial positions of the Galilean 
moons are indicated at the bottom of the figure. 


can be inferred more easily from observations than age. Hence, 
planetary accretion rates allow us to consider planet (and moon) 
formation at comparable stages, independent of the initial con¬ 
ditions. 

Variations of the stellar distance will affect both the size of 
the circumplanetary accretion disk (Machida et al. 2008) and the 
stellar heating term in our model (Makalkin & Dorofeeva 2014). 
Disks around close-in planets are usually smaller (though not 
necessarily less massive) than disks around planets in wide stel¬ 
lar orbits owing to their smaller Hill spheres and their lower av¬ 
erage specific angular momentum. Moreover, substantial stellar 
illumination will prevent these disks from having HjO ice lines 
in the vicinity of the star. The temperature of the circumstellar 
accretion disk, in which the planet is embedded, is calculated 
under the assumption that the disk is transparent to stellar irradi¬ 
ation (Hayashi 1981) and that the stellar luminosity equals that 
of the contemporary Sun. In this model, the circumstellar H 2 O 
ice line is located at 2.7 AU from the star. We will revisit the faint 
young Sun (Sagan & Mullen 1972) as well as other parameter- 
izations of the circumstellar disk in an upcoming paper (Heller 
et al. 2015, in prep.). 

In Figs. 2 and 3 we show an application of our disk model 
to a Jupiter-mass and a 12Mj up planet at 5.2 AU around a Sun¬ 
like star, respectively. All panels present the disk surface tem¬ 
peratures ( T s , solid gray lines) and disk midplane temperatures 
(T m , solid black lines) as well as the contributions to T s from 
viscous heating ( T v , black crosses), planetary illumination 17),, 
solid black lines), accretion onto the disk ( 7' a , gray dashed lines), 
and heating from the circumstellar accretion disk, or “nebula” 
(7’ n , gray dotted lines). Any of these contributions to the disk 
surface temperature is computed assuming that all other contri¬ 
butions are zero. In other words, T v , T p , T a , and 7' n depict the 
temperature of the disk photosphere in the hypothetical case that 
viscous heating, planetary illumination, accretional heating onto 
the disk, or the stellar illumination were the single energy source, 
respectively. 3 The different slopes of these curves, in particu¬ 
lar of T v and T v , lead to the appearance of heat transitions (not 
shown), which traverse the disk on relatively short time scales 


3 As an example, T v is calculated setting F acc = F p = T neb = 0 in 
Eq. (13) of Heller & Pudritz (2014). 


(Heller & Pudritz 2014). Also note that in all the simulations 
shown, the disk midplane is warmer than the disk surface. Only 
when accretion drops to about 10 Mq w Myr - in panels (b) do 
T s and T m become comparable throughout the disk, both around 
the Jupiter- and the super-Jupiter-mass planet. 

Panels (a) in Figs. 2 and 3 are chosen at the time at which 
the H 2 O ice line first appears at the outer edge of the disk, 
while panels (b) illustrate the temperature distribution at the time 
when M = 10 Mean, with Moan as Ganymede’s mass. During 
this epoch, the H 2 O ice line around the Jupiter-mass planet is 
safely between the orbital radii of Europa and Ganymede (see 
labeled arrow in Fig. 2b), which suggests that this corresponds 
to the shutdown phase of moon formation. The colored symbols 
in each panel depict the rocky (orange) or icy (blue) composi¬ 
tion of the Galilean moons. Symbol sizes scale with the actual 
moon radii. In these simulations, the Planck mean opacity ( Kp ) 
has been fixed at a fiducial value of 10 2 m 2 kg -1 throughout the 
disk and the disk absorptivity is assumed to be k s - 0.2. 

A comparison of Figs. 2(a) and 3(a) shows that the cen¬ 
trifugal radius of the accretion disk around the Jovian planet 
(r c ~ 27 Fj up ) is substantially smaller than the disk radius around 
the super-Jovian test planet (r c ~ 807?j up ) - note the different 
distance ranges shown in the figures! Moreover, the midplanes 
and surfaces in the inner regions of the super-Jovian accretion 
disks are substantially hotter at any given distance around the 
super-Jupiter. Note also that planetary illumination is the main 
energy source around the Jupiter-mass planet in Figs. 2(a) and 
(b), while viscous heating plays an important role during the ap¬ 
pearance of the H 2 O ice line around the super-Jovian planet in 
Fig. 3(a). In the final stages of accretion onto the 12 Mj up planet, 
shown in Fig. 3(b), viscous heating and planetary illumination 
are comparable throughout the disk. 

In Fig. 4 we present the radial distributions of the gas sur¬ 
face densities around our Jovian and super-Jovian test planets. 
While solid lines refer to panels (a) in Figs. 2 and 3, dashed 
lines refer to panels (b), respectively. In particular, solid lines 
refer to that instant in time when the H 2 O ice lines first appear 
at the outer edges of the accretion disks around those two plan¬ 
ets, whereas the dashed lines depict the accretion phase when 
M - 10 Mean Myr -1 , which is the phase when the H 2 O ice line 
around our Jupiter-mass test planet is in good agreement with the 
compositional gradient observed in the Galilean system (Heller 
& Pudritz 2014). At any given planetary distance in these partic¬ 
ular states, the super-Jovian planet has a gas surface density that 
is higher by about a factor of five compared to the Jupiter-mass 
planet. Moreover, we find that the H 2 O ice lines in both cases 
need about 10 s yr to move radially from the outer disk edges to 
their final positions (see labels). Our values for X g are similar to 
those presented by Makalkin & Dorofeeva (2014) in their Fig. 4. 

Similar to the procedure applied in Heller & Pudritz (2014), 
we performed a suite of simulations for planets with masses be¬ 
tween 1 and 12 Mj up at 5.2 AU from a Sun-like star, where M s h u t 
and Kp were randomly drawn from a Gaussian probability dis¬ 
tribution. For log 10 (Kp/[m 2 kg -1 ]) we took a mean value of -2 
with a standard variation of 1 (see Bell et al. 1997, and refer¬ 
ences therein), and for log 10 (M s h u t/[A7Gan Myr -1 ]) we assumed a 
mean value of 1 with a standard variation of 1, which nicely re¬ 
produced the compositional H 2 O gradient in the Galilean moons 
(Heller & Pudritz 2014). As an extension of our previous simula¬ 
tions, we here consider various disk absorptivities (k s = 0.2 and 
0.4) and examine the total mass of solids in the accretion disks, 
both as a function of time and as a function of stellar distance. 
Above all, we study the disappearance of the circumplanetary 
H 2 O ice line in the vicinity of the star. 
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Fig. 5. Distance of the circumplanetary H 2 O ice lines as a function of planetary mass. A distance of 5.2 AU to a Sun-like star is assumed. The 
shaded area indicates the lcr scatter of our simulations based on the posterior distribution of the disk Planck mean opacity (k p ) and the shutdown 
accretion rate for moon formation (M s hut). The labeled circles at 1 Mj llp denote the orbital positions of Jupiter's moons Io, Europa, Ganymede, and 
Callisto. Orange indicates rocky composition, blue represents TEO-rich composition. Circle sizes scale with the moons’ radii, (a): Disk reflectivity 
( k s ) is set to 0.2. (b): Same model parameterization but now with disk reflectivity k s = 0.4. 




Fig. 6. Similar to Fig. 5, but now in units of fractional disk radius (see ordinate), (a): Disk reflectivity (k s ) is set to 0.2. (b): Same model parameter¬ 
ization but now with disk reflectivity k s = 0.4. The negative slope in both panels indicates that the most massive giant planets have larger fractions 
of their disks fed with water ice and, consequently, relatively more space and material to form big moons. 


3. H 2 0 ice lines around planets of different masses 

Figures 2(a) and (b) show that the temperature distribution in the 
late accretion disks around Jupiter-mass planets is determined 
mostly by the planetary illumination rather than viscous heat¬ 
ing. This indicates a key difference between moon formation in 
circumplanetary accretion disks and planet formation in circum- 
stellar disks, where the positions of ice lines have been shown 
to depend mostly on viscous heating rather than stellar illumina¬ 
tion (Min et al. 2011; Hasegawa & Pudritz 2011). In this con¬ 
text, Figs. 3(a) and (b) depict an interesting intermediate case 
in terms of the mass of the central object, where both viscous 
heating and illumination around a 12Mj up planet show compa¬ 
rable contributions towards the final stages of accretion in panel 
(b). Also note that in Fig. 3(b) the temperature distribution in 
the disk outskirts is determined by the background temperature 
provided by stellar illumination. This suggests that the extended 
disks around super-Jovian exoplanets in wide circumstellar or¬ 
bits (beyond about 10 AU) might also feature CO and other ice 
lines due to the even weaker stellar illumination. 

Figures 5(a) and (b) display the radial positions of the cir¬ 
cumplanetary H 2 O ice lines around a range of planets with 
masses between 1 and 12Mj up in the final stages of accretion. 


where we randomized Kp and M s h u t as described above, and the 
stellar distance is 5.2 AU for all planets considered. Panel (a) 
shows the same data as Fig. 4 in Heller & Pudritz (2014), while 
panel (b) assumes k s - 0.4. The larger absorptivity in panel (b) 
pushes the H 2 O ice lines slightly away from the accreting plan¬ 
ets compared to panel (a) for all planets except for the 12 Mj up 
planet, which we ascribe to an insignificant statistical fluctua¬ 
tion. As a key result of our new simulations for k s = 0.4, and in 
agreement with our previous study, the H 2 O ice line is between 
about 15 and 30 Rj ap for all super-Jovian planets and almost in¬ 
dependent of M p . The effect of changing k s by a factor of two is 
moderate, pushing the H 2 O ice line outward by only a few Rj up 
on average. 

Figure 6 presents a different visualization of these results, 
again with panel (a) assuming k s = 0.2 and panel (b) assuming 
k s = 0.4. Now the ordinate gives the circumplanetary distance in 
units of the fractional disk radius, and so the H 2 O ice lines are 
located closer to the planet with increasing M p . This is a conse¬ 
quence of the larger disk sizes of the more massive planets. In 
other words, larger giant planets have larger fractions of the ac¬ 
cretion disks beyond the H 2 O ice lines. Naturally, this means that 
these super-Jovians should form the most massive, icy moons. 
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Fig. 7. Evolution of the total mass of solids (M sld ) in circumplane- 
tary accretion disks of several test planets at 5.2 AU from a Sun-like 
star. Colored dots along the mass evolution tracks indicate accretion 
rates (M). If moon formation around all of these planets shuts down at 
comparable disk accretion rates (A/ shut ), then M ski (M shllt ) oc M p (see also 
Fig. 5 in Heller & Pudritz 2014). Compared to a Jupiter-mass planet, a 
12 Af Jup planet would then have 12 times the amount of solids available 
for moon formation. 

In Fig. 7, we therefore analyze the evolution of M s y (in units 
of Moan) around these planets as a function of time. Indicated 
values for M along these tracks help to visualize M s \& as a func¬ 
tion of accretion rates (see crosses, squares, and circles), because 
time on the abscissa delivers an incomplete picture of the mass 
evolution due to the different timescales of these disks. We find 
that disks around lower-mass super-Jovians contain less solid 
mass at any given accretion rate than disks around higher-mass 
super-Jovians. Solids also occur earlier in time as the appearance 
of the disk itself is regulated by the shrinking of the initially 
very large planet in our models: lower-mass giant planets con¬ 
tract earlier (Mordasini 2013). Assuming that moon formation 
shuts down at similar accretion rates around any of the simulated 
planets, the example accretion rates (100, 10, 1 McanMyU 1 ) in¬ 
dicate an increase of the mass of solids as a function of M p . For 
a given M s h u t, this scaling is M s y oc M p (see also Fig. 5 in Heller 
& Pudritz 2014), which is in agreement with the scaling relation 
for the total moon mass around giant planets found by Canup & 
Ward (2006). 

4. H 2 0 ice lines and total mass of solids around 
planets at various stellar distances 

Having examined the sensitivity of our results to M p and to 
the disk’s radiative properties, we now turn to the question of 
what exomoon systems are like around super-Jovian exoplanets 
at very different locations in their disks than our own Jupiter at 
5.2 AU. Here, we shall encounter some significant surprises. 

Figure 8(a) shows the final circumplanetary distance of 
the H 2 0 ice line around a Jupiter-mass planet between 2 and 
20 AU from the star, 4 assuming a shutdown accretion rate of 
lOOMcanMyr 1 . Different styles of the blue lines correspond 
to different disk opacities (see legend), while the solid black 

4 We also simulated planets as close as 0.2 AU to the star, but their ac¬ 
cretion disks are nominally smaller than the planetary radius, indicating 
a departure of our model from reality. Anyways, since the hypothetical 
disks around these planets do not harbor H 2 0 ice lines and moon for¬ 
mation in the stellar vicinity is hard in the first place (Barnes & O’Brien 
2002; Namouni 2010), we limit Fig. 10 to 2 AU. 


line indicates the radius of the circumplanetary accretion disk, 
following Machida et al. (2008). Circles indicate the radial dis¬ 
tances of the Galilean moons around Jupiter, at 5.2 AU from the 
star. In this set of simulations, the H 2 0 ice lines always ends 
up between the orbits of Ganymede and Callisto, which is not 
in agreement with the observed H 2 0 compositional gradient in 
the Galilean system. As stellar illumination decreases at larger 
distances while all other heating terms are constant for the given 
accretion rate, the ice lines move towards the planet at greater 
distances. 

Most importantly, however, we find that Jovian planets closer 
than about 4.8 AU do not have an H 2 0 ice line in the first place. 
Hence, if the large population of Jupiter-mass planets around 
1 AU (see Fig. 1) formed in situ and without substantial in¬ 
ward migration from the outer regions, then these giant plan¬ 
ets should not have had the capacity to form giant, icy moons, 
that is, scaled-up versions of Ganymede or Callisto. These giant 
moons with masses up to that of Mars (suggested by Canup & 
Ward 2006; Heller et al. 2014; Heller & Pudritz 2014), may only 
be present if they have completed their own water-rich forma¬ 
tion beyond about 4.8 AU from their star, before they migrated 
to their current circumstellar orbits together with their host plan¬ 
ets. 

In Fig. 8(b), the planetary accretion rate has dropped 
by a factor of ten, and the ice lines have shifted. For 

kp — 1 (J- 2 m 2 kg 1 and Kp - 10 1 nr kg 1 they moved towards 

the planet. But for Kp = 10~ 3 m 2 kg 1 they moved outward. 

The former two rates actually place the H 2 0 ice line around 
Jupiter at almost exactly the orbit of Ganymede, which is in 
better agreement with observations. In these simulations, Jovian 
planets closer than about 4.5 AU do not have a circumplanetary 
H 2 0 ice line. 

In Fig. 9 we vary the stellar distance of a 12Mj up planet. 
First, note that the disk (black solid line) is larger at any given 
stellar separation than in Fig. 8. Second, note that the H 2 0 ice 
line for a given M s h u t in Figs. 9(a) and (b), respectively, is fur¬ 
ther out than in the former case of a Jupiter-mass planet. This is 
due to both increased viscous heating and illumination from the 
planet for the super-Jovian object. Nevertheless, the much larger 
disk radius overcompensates for this effect and, hence, all these 
simulations suggests that accretion disks around the most mas¬ 
sive planets can have H 2 0 ice lines if the planet is not closer than 
about 3.8 AU (panel a) to 3.1 AU (panel b) to a Sun-like star. The 
latter value refers to a shutdown accretion rate of 10 Mean MyU 1 , 
which is in good agreement with the water ice distribution in the 
Galilean moon system. Assuming that moon formation stops at 
comparable accretion rates around super-Jovian planets, we con¬ 
sider a critical stellar distance of about 3 AU a plausible estimate 
for the critical stellar distance of a 12Mj up accreting planet to 
show a circumplanetary H 2 0 ice line. 

Figure 10 presents the total instantaneous mass of solids in 
circumplanetary accretion disk as a function of stellar distance 
for two different shutdown accretion rates and three different 
disk Planck opacities. The calculation of M s y follows Eq. (17) 
in Heller & Pudritz (2014), that is, we integrate the surface 
density of solids between the inner disk truncation radius and 
the outer centrifugal radius. Panel (a) for a Jupiter-mass planet 
demonstrates that indeed the amount of solids in its late accretion 
disk is negligible within about 4.5 AU from a Sun-like star. This 
gives us crucial insights into Jupiter’s migration history within 
the Grand Tack framework (Walsh et al. 2011), which we will 
present in a forthcoming paper (Heller et al. 2015, in prep.). 
Moreover, both panels show that the effect of different disk opac¬ 
ities on M s id is small for a given M p and M s h u t- 
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Fig. 8. Distances of the circumplanetary H 2 O ice lines around a Jupiter-like planet (ordinate) as a function of distance to a Sun-like star (abscissa) 
at the time of moon formation shutdown. Each panel assumes a different shutdown formation rate (see panel titles). Three Planck opacities through 
the circumplanetary disk are tested in each panel (in units of m 2 kg -1 , see panel legends). The circumplanetary orbits of the Galilean satellites are 
represented by symbols as in Fig. 2. The black solid line shows the disk’s centrifugal radius (Machida et al. 2008). (a) At M shut = 100 M Gm /Myr the 
H 2 O ice lines at » 5 AU are about 5 Rj up beyond Ganymede’s current orbit, (b) At M S h U t = 10 Moan/Myr values of 1(T 2 m 2 kg~* < k < 10 _1 m 2 kg~* 
place the HiO ice line slightly inside the current orbit of Ganymede and thereby seem most plausible. 


12 M Jup M shut =100M Gan /Myr 12M Jup M shut = 10M Gim /Myr 




Fig. 9. Same as Fig. 8, but now for a 12 Mj up mass planet. 


Intriguingly, for any given disk parameterization (see legend) 
a circumjovian disk at 5.2 AU from a Sun-like star seems to har¬ 
bor a relatively small amount of solids compared to disks around 
planets at larger stellar distances. The red dashed line indicates 
a best fit exponential function to the simulations of our fiducial 
disk with M s hut = 10 Myr 1 and Kp = 10 _2 m 2 kg 1 beyond 
5.2 AU, as an example. We choose this disk parameterization as 
it yields the best agreement with the radial location of the icy 
Galilean satellites (Heller & Pudritz 2014). It scales as 

M sld = 10~ 5 - 3 M Gan x (^y) 1 ' 16 , (1) 

where r* is the stellar distance. Assuming that moon formation 
stops at similar mass accretion rates around any super-Jovian 
planet and taking into account the previously known scaling of 
the total moon masses ( Mj ) with M p , we deduce a more general 
estimate for the total moon mass: 




where Mgm = 2.65 Mean is the total mass of the Galilean moons. 
As an example, a 10 Mj up planet forming at 5.2 AU around a Sun¬ 
like star should have a moon system with a total mass of about 
IOMqm or 6 times the mass of Mars. At a Saturn-like stellar 
distance of 9.6 AU, Mr would be doubled. Simulations by Heller 
et al. (2014) show that this mass will be distributed over three to 
six moons in about 90 % of the cases. Hence, if the most massive 
super-Jovian planets formed moon systems before they migrated 
to about 1 AU, where we observe them today (see Fig. 1), then 
Mars-mass moons in the stellar habitable zones might be very 
abundant. 


5. Shutdown accretion rates and loss of moons 

Canup & Ward (2002) argued that accretion rates of 
2 x 10~ 7 Mj up yr -1 (or about 2.6 x 10 4 Mean MyU 1 ) best repro¬ 
duce the disk conditions in which the Galilean system formed. 
Based on the condition that the H 2 O ice line needs to be be¬ 
tween the orbits of Europa and Ganymede at the final stages 
of accretion, our calculations predict a shutdown accretion rate 
that is considerably lower, closer to 10 Mean Myr 1 (see Fig. 8). 
The difference in these results is mainly owed to two facts. 
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Fig. 10. Total instantaneous mass of solids in circumplanetary accretion disks in units of Ganymede masses and as a function of stellar distance. 
Labels indicate two different values for M s hut, and different line types refer to different disk opacities (see legend). Black dashed lines in both 
panels indicate our fiducial reference disk, (a): The 1 M Jllp test planet starts to have increasingly large amounts of solids only beyond 4.5 AU, 
depending on the accretion rate and disk opacity. The red dashed line indicates our fit as per Equation (1). (b): The 12 Mj up test planet starts to 
contain substantial amounts of solids as close as about 3.1 AU to the star, owed to its larger disk. 


First, Canup & Ward (2002) only considered viscous heating and 
planetary illumination. Our additional heating terms (accretion 
onto the disk and stellar illumination) contribute additional heat, 
which imply smaller accretion rates to let the FLO ice lines move 
close enough to the Jupiter-like planet. Second, the parameteri¬ 
zation of planetary illumination in the Canup & Ward (2002) 
model is different from ours. While Canup & Ward (2002) as¬ 
sume an r !/4 dependence of the midplane temperature from the 
planet (r being the planetary radial distance), we do not apply 
any predescribed r-dependence. In particular, T m (r) cannot be 
described properly by a simple polynomial due to the different 
slopes of the various heat sources as a function of planetary dis¬ 
tance (see the black solid lines in Figs. 2 and 3). 

While our estimates of M s h u t are about three orders of mag¬ 
nitude lower than the values proposed by Canup & Ward (2002), 
they are also one to two orders of magnitude lower than the val¬ 
ues suggested by Alibert et al. (2005). Their moon formation 
model is similar to the so-called “gas-starved” disk model ap¬ 
plied by Canup & Ward (2002). Accretion rates in their model 
were not derived from planet formation simulations (as in our 
case) but calculated using an analytical fit to previous simula¬ 
tions. Again, the additional energy terms in our model are one 
reason for the lower final accretion rates that we require in order 
to have the FLO ice line interior to the orbit of Ganymede. 

Makalkin & Dorofeeva (2014) also studied Jupiter’s ac¬ 
cretion rates and their compatibility with the Galilean moon 
system. They found that values between 1 0 -9 Mj up yr _I and 
10 _6 M Jup yU 1 (about lOMcanMyU 1 to 10 4 Mc a nMyU 1 ) satisfy 
these constraints. Obviously, our results are at the lower end 
of this range, while accretion rates up to three orders of mag¬ 
nitude higher should be reasonable according to Makalkin & 
Dorofeeva (2014). Yet, these authors also claimed that plane¬ 
tary illumination be negligible in the final states of accretion, 
which is not supported by our findings (see Fig. 2). We ascribe 
this discrepancy to the fact that they estimated the planetary lu¬ 
minosity (L p ) analytically, while in our model L p , R v , M p , and 
M are coupled in a sophisticated planet evolution model (Mor- 
dasini 2013). In the Makalkin & Dorofeeva (2014) simulations, 
10~ 7 Lq < L p < 1CL 4 L 0 (depending on assumed values for 
R p , M p , and M), whereas in our case L p remains close to 10 4 L 0 
once a circumjovian FLO ice line forms, while M and R p evolve 


rapidly due to the planet’s gap opening in the circumsolar disk 
(see Fig. 1 in Heller & Pudritz 2014). 

We also calculate the type I migration time scales (ti) of 
potential moons that form in our circumplanetary disks. Using 
Eq. (1) from Canup & Ward (2006) and assuming a Ganymede- 
mass moon, we find that 0.1 Myr < n < 100 Myr in the 

final stages of accretion, with the shortest time scales referring 
to close-in moons and early stages of accretion when the gas sur¬ 
face density is still high. Hence, since the remaining disk lifetime 
(ss 10 6 yr) is comparable or smaller than the type I migration 
time scale, migration traps might be needed to stop the moons 
from falling into the planet. As the gas surface density is de¬ 
creasing by about an order of magnitude per 10 5 yr (see Fig. 1(c) 
in Heller & Pudritz 2014) and since ti is inversely proportional 
to Eg, type I migration slows down substantially even if the pro¬ 
tosatellites still grow. On the other hand, if a moon grows large 
enough to open up a gap in the circumplanetary accretion disk, 
then type II (inwards) migration might kick in, reinforcing the 
need for moon migration traps. All these issues call for a de¬ 
tailed study of moon migration under the effects of ice lines and 
other traps. 

Another issue that could cause moon loss is given by their 
possible tidal migration. If a planet rotates very slowly, its coro¬ 
tation radius will be very wide and any moons would be forced 
to tidally migrate inwards (ignoring mean motion resonances for 
the time being). But if the planet rotates quickly and is not sub¬ 
ject to tides raised by the star, such as Jupiter, then moons usually 
migrate outwards due to the planetary tides and at some point 
they might become gravitationally unbound. Barnes & O’Brien 
(2002, see their Fig. 2) showed that Mars-mass moons around gi¬ 
ant planets do not fall into the planet and also remain bound for 
at least 4.6 Gyr if the planet is at least 0.17 AU away from a Sun¬ 
like star. Of course, details depend on the exact initial orbit of the 
moon and on the tidal parameterization of the system, but as we 
consider planets at several AU from the star, we conclude that 
loss of moons due to tidal inward or outward migration is not an 
issue. Yet, it might have an effect on the orbital distances where 
we can expect those giant moons to be found, so additional tidal 
studies will be helpful. 
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6. Discussion 

If the abundant population of super-Jovian planets at about 1 AU 
and closer to Sun-like stars formed in situ, then our results sug¬ 
gest that these planets could not form massive, super-Ganymede- 
style moons in their accretion disks. These disks would have 
been too small to feature hkO ice lines and therefore the growth 
of icy satellites. The accretion disks around these planets might 
still have formed massive, rocky moons, similar in composition 
to Io or Europa, which would likely be in close orbits (thereby 
raising the issue of tidal evolution, Barnes & O’Brien 2002; Cas¬ 
sidy et al. 2009; Porter & Grundy 2011; Heller & Barnes 2013; 
Heller et al. 2014), since circumplanetary accretion disks at 1 AU 
are relatively small. If these large, close-orbit rocky moons exist, 
they also might be subject to substantial tidal heating. Alterna¬ 
tively, super-Jovians might have captured moons, e.g. via tidal 
disruption of binary systems during close encounters (Agnor & 
Hamilton 2006; Williams 2013), so there might exist indepen¬ 
dent formation channels for giant, possibly water-rich moons at 
1 AU. 

If, however, these super-Jovian planets formed beyond 3 to 
4.5 AU and then migrated to their current locations, then they 
could be commonly orbited by Mars-mass moons with up to 
50 % of water, similar in composition to Ganymede and Callisto. 
Hence, the future detection or non-detection of such moons will 
help to constrain rather strongly the migration history of their 
host planets. What is more, Mars-mass ocean moons at about 
1 AU from Sun-like stars might be abundant extrasolar habitats 
(Williams et al. 1997; Heller & Barnes 2013; Heller et al. 2014), 

Our results raise interesting questions about the formation 
of giant planets with satellite systems in the solar system and 
beyond. The field of moon formation around super-Jovian exo¬ 
planets is a new research area, and so many basic questions still 
need to be answered. 

(1) The “Grand Tack”. In the “Grand Tack” scenario 
(Walsh et al. 2011), the fully accreted Jupiter migrated as close 
as 1.5 AU to the Sun during the first few million years of the 
solar system, then got caught in a mean-motion orbital reso¬ 
nance with Saturn and then moved outward to about 5 AU. Our 
results suggest that the icy moons Ganymede and Callisto can 
hardly have formed during the several 10 5 yr Jupiter spent inside 
4.5 AU to the Sun. If they formed before Jupiter’s tack, could 
their motion through the inner solar system be recorded in these 
moons today? Alternatively, if Callisto required 10 5 - 10 6 yr to 
form (Canup & Ward 2002; Mosqueira & Estrada 2003) and 
assuming that Jupiter’s accretion disk was intact until after the 
tack, one might suppose that Ganymede and Callisto might have 
formed thereafter. But then how did Jupiter’s accretion disk re¬ 
acquire the large amounts of H 2 O that would then be incorpo¬ 
rated into Callisto after all water had been sublimated during the 
tack? We will address these issues in a companion paper (Heller 
et al. 2015). 

(2) Migrating planets. Future work will need to include the 
actual migration of the host planets, which we neglected in this 
paper. Namouni (2010) studied the orbital stability of hypothet¬ 
ical moons about migrating giant planets, but the formation of 
these moons was not considered. Yet, the timing of the accretion 
evolution, the gap opening, the movement within the circumstel- 
lar disk (and thereby the varying effect of stellar heating), and 
the shutdown of moon formation will be crucial to fully assess 
the possibility of large exomoons at about 1 AU. These simula¬ 
tions should be feasible within the framework of our model, but 
the precomputed planet evolution tracks would need to consider 
planet migration. Ultimately, magneto-hydrodynamical simula¬ 


tions of the circumplanetary accretion disks around migrating 
super-Jovian planets might draw a full picture. 

(3) Directly imaged planets. Upcoming ground-based ex¬ 
tremely large telescopes such as the E-ELT and the Thirty Me¬ 
ter Telescope, as well as the James Webb Space Telescope have 
the potential to discover large moons transiting directly imaged 
planets in the infrared (Peters & Turner 2013; Heller & Albrecht 
2014). Exomoon hunters aiming at these young giant planets be¬ 
yond typically 10 AU from the star will need to know how moon 
formation takes place in these possibly very extended circum¬ 
planetary accretion disks, under negligible stellar heating, and in 
the low-density regions of the circumstellar accretion disk. 

(4) Ice line traps. If circumplanetary HjO ice lines can act 
as moon migration traps and if solids make up a substantial part 
of the final masses accreted by the planet, then the accretion rates 
onto giant planets computed under the neglect of moons might 
be incorrect in the final stages of accretion. The potential of the 
H 2 O ice line to act as a moon migration trap is new (Heller & Pu- 
dritz 2014) and needs to be tested. It will therefore be necessary 
to compute the torques acting on the accreting moons within the 
circumplanetary accretion disk as well as the possible gap open¬ 
ing in the circumplanetary disk by large moons, which might 
trigger type II migration (Canup & Ward 2002, 2006). 


7. Conclusions 

Planetary illumination is the dominant energy source in the late- 
stage accretion disks around Jupiter-mass planets at 5.2 AU from 
their Sun-like host stars (Fig. 2), while viscous heating can be 
comparable in the final stages of accretion around the most mas¬ 
sive planets (Fig. 3). 

At the time of moon formation shutdown, the H 2 O ice line 
in accretion disks around super-Jovian planets at 5.2 AU from 
Sun-like host stars is between roughly 15 and 30/?j up . This dis¬ 
tance range is almost independent of the final planetary mass and 
weakly dependent on the disk’s absorption properties (Fig. 5). 
With more massive planets having more extended accretion 
disks, this means that more massive planets have larger fractions 
(up to 70 %) of their disks beyond the circumplanetary H 2 O ice 
line (compared to about 25 % around Jupiter, see Fig. 6). 

Jupiter-mass planets forming closer than about 4.5 AU to a 
Sun-like star do not have a circumplanetary H 2 O ice line (Fig. 8), 
depending on the opacity details of the circumstellar disk. A de¬ 
tailed application of this aspect to the formation of the Galilean 
satellites might help constraining the initial conditions of the 
Grand Tack paradigm. Due to their larger disks, the most mas¬ 
sive super-Jovian planets can host an H 2 O ice line as close as 
about 3 AU to Sun-like stars (Fig. 9). With the circumstellar H 2 O 
ice line at about 2.7 AU in our model of an optically thin cir¬ 
cumstellar disk (Hayashi 1981), the relatively small accretion 
disks around Jupiter-mass planets at several AU from Sun-like 
stars thus substantially constrain the formation of icy moons. 
The extended disks around super-Jovian planets, on the other 
hand, might still form icy moons even if the planet is close to the 
circumstellar H 2 O ice line. 

We find an approximation for the total mass available for 
moon formation (Mr), which is a function of both M p and r* (see 
Eq. 2). The linear dependence of Mp °c M p has been known 
before, but the dependence on stellar distance (Mp oc 
for r* > 5.2 AU) is new. It is based on our finding that an 
accretion rate M s hut ~ 10 Mean/ 10 6 yr (Fig. 10) yields the best 
results for the position of the circumjovian H 2 O ice line at the 
shutdown of the formation of the Galilean moons (Fig. 2b); and 
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it assumes that this shutdown accretion rate is similar around all 
super-Jovian planets. 

Our results suggest that the observed large population of 
super-Jovian planets at about 1 AU to Sun-like stars should not 
be orbited by water-rich moons if the planets formed in-situ. 
However, in the more plausible case that these planets migrated 
to their current orbits from beyond about 3 to 4.5 AU, they 
should be orbited by large, Mars-sized moons with astrobiolog- 
ical potential. As a result, future detections or non-detections of 
exomoons around giant planets can help to distinguish between 
the two scenarios because they are tracers of their host planets’ 
migration histories. 
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